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Abstract: We report on the status of supersymmetric seesaw models in the light of 
recent experimental results on fi — > e + 7, #13 and the light Higgs mass at the LHC. 
SO(10)-like relations are assumed for neutrino Dirac Yukawa couplings and two cases of 
mixing, one large, PMNS-like, and another small, CKM-like, are considered. It is shown 
that for the large mixing case, only a small range of parameter space with moderate tan f3 is 
still allowed. This remaining region can be ruled out by an order of magnitude improvement 
in the current limit on BR(/x — > e + 7). We also explore a model with non- universal Higgs 
mass boundary conditions at the high scale. It is shown that the renormalization group 
induced flavor violating slepton mass terms are highly sensitive to the Higgs boundary 
conditions. Depending on the choice of the parameters, they can either lead to strong 
enhancements or cancellations within the flavor violating terms. Such cancellations might 
relax the severe constraints imposed by lepton flavor violation compared to mSUGRA. 
Nevertheless for a large region of parameter space the predicted rates lie within the reach 
of future experiments once the light Higgs mass constraint is imposed. We also update 
the potential of the ongoing and future experimental searches for lepton flavor violation in 
constraining the supersymmetric parameter space. 
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1 Introduction 

Current times are unprecedented in terms of experimental activity in high energy physics. 
There have been several results of very high impact in the recent times. The following 
three are the most relevant for the purpose of our discussion. 

• Firstly, LHC experiments have reported the discovery of a new boson with mass 
of about 125 GeV [1, 2], compatible with the Standard Model (SM) Higgs boson. 
In minimal supersymmetric standard model (MSSM), this would imply a mass of 
this order and SM-like couplings for the lightest CP-even Higgs boson [3-6]. In our 
analysis, we will take the 2a mass range obtained in [3] 

124.5 GeV < m h < 126.5 GeV, (1.1) 

where stands for the mass of the lightest neutral Higgs. In addition, LHC has 
also improved the limits on the spectrum of low-energy supersymmetry (SUSY) [7] . 

• The limit on the lepton flavor violating (LFV) decay ji — > e + 7 has improved by one 
order of magnitude [8]. The current limit is 

BRGu -> e + 7) < 2.4 x 1(T 12 (90% CL). (1.2) 

• Finally, the so-far unknown neutrino mixing angle, #13 has been experimentally de- 
termined [9, 10]. The Daya Bay and RENO experiments have measured #13 with 
very good accuracy: 

sin 2 26»i3 = 0.092 ± 0.016(stat.) ± 0.005(syst.) Daya Bay [9] 

sin 2 20i3 = 0.113 ± 0.013(stat.) ± 0.019(syst.) RENO [10] (1.3) 
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The implications of the measurement of the light Higgs mass for various supersymmet- 
ric models have been studied in detail by various authors [11-23]. Supersymmetric seesaw 
models have not been explored so far in the light of these results 1 . Furthermore, flavor 
violation is expected from SUSY seesaw models where the last two experimental results 
would play a crucial role in constraining the parameter space. In the present work, we 
consider the implications of all the three experimental results on a class of SUSY seesaw 
models inspired by SO(10) GUTs. 

SO (10) GUT Models typically relate up-quark Yukawa matrices with the Dirac neu- 
trino Yukawa couplings of the Type I seesaw mechanism. In fact, at least one of the 
neutrino Yukawa couplings is expected to be as large as the top Yukawa coupling, as a 
consequence of the undelying SO (10) symmetry [25]. Concerning the mixing structure of 
the neutrino Yukawa matrix, two extreme cases can be motivated by simple SO(10) models: 
large PMNS-like mixing or small CKM-like mixing. Lepton flavor violation in this class 
of models has been studied in Refs. [25-27]. The present work can be considered as an 
update of these studies. 

The previous works considered mSUGRA/CMSSM-type boundary conditions for the 
SUSY-breaking soft terms. While universality of the soft terms at the GUT scale is re- 
quired for the fields belonging to the same 16 representation of SO(10), the full universal 
boundary conditions of the CMSSM type are too strong a condition. For example, there 
is no fundamental reason why the Higgs doublets (that typically are in 10-dimensional 
representations) and the sfermion soft masses should be degenerate. Thus, strict univer- 
sality between matter and Higgs fields can be relaxed [28]. These models are typically 
dubbed Non- Universal Higgs Mass models (NUHM). The relaxation of the universality has 
an important impact on the seesaw generated flavor violating entries of the slepton mass 
matrix. The magnitude of the RG generated flavor violating entries can either increase 
or decrease at the leading order due to the interplay between the matter and Higgs mass 
terms at the GUT scale. Cancellations in the flavor violating entry can indeed relax the 
LFV constraints on the SUSY parameter space. However, as we are going to see, the Higgs 
mass and LFV constraints are such that we have similar conclusions as in the mSUGRA, 
for moderate/large values of tan/3: an improvement of one order of magnitude in the 
BR(/x — t- e + 7) bound is sufficient to rule out significant amount of the parameter space. 
In the present work, we compare and contrast the constraints on SUSY seesaw parameter 
space with CMSSM/mSUGRA boundary conditions and NUHM boundary conditions. 

We find that in the PMNS case, tan/3 is restricted between 4 to 20. Lower values of 
tan (3 are disfavored by the light Higgs mass constraints whereas higher values are strongly 
constrained by the present limit on BR(^ — > e + 7) whose rates can be very large, as a 
consequence of the sizable observed value of #13. Furthermore, the one order of magnitude 
improvement expected in the near future from the MEG experiment would rule out most 
of the SUSY parameter space accessible at the LHC. This situation is somewhat relaxed 
in the NUHM model where we considered the masses of both the Higgs doublets to be the 
same (the so-called NUHM1 model). We have further updated the future prospects for 
both the CKM and PMNS cases in mSUGRA/CMSSM as well as in NUHM1. 

The paper is organized as follows. In section 2, we recap the SUSY seesaw and discuss 
the generation of flavor violation in mSUGRA and NUHM1. In section 3, we discuss the 

1 As this work was finished and being prepared for submission, the following paper Ref. [24] appeared on 
arXiv. 
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details of our numerical analysis. In section 4, we present our results. We conclude with a 
summary and outlook in section 5. Finally, in appendix A we describe the proposed future 
experiments and their expected sensitivity. 



2 Seesaw in mSUGRA and NUHM 

The phenomenology of SUSY Type I seesaw mechanism with universal boundary conditions 
(mSUGRA/CMSSM) has been studied in many papers (see [29, 30] for a set of recent 
works). Here we review some essential features related to flavor violation for completeness 
and to do a comparison with the case of non-universal Higgs masses. To set the notation, 
the Type I seesaw mechanism is characterized by a superpotential containing the following 
terms 

W D Y e Le c H d + Y v Lv c H u + ~M R v c v c (2.1) 

where L (e c ) stands for the leptonic doublets (singlets) and v c are the right-handed (RH) 
neutrino superfields (with the generation indices not explictely written). Y e and Y„ are 
the electron and neutrino (Dirac) Yukawa matrices. 

In models like CMSSM/mSUGRA, the soft terms are assumed to be universal at the 
Grand Unification (GUT) scale, Mgut ~ 2 x 10 16 GeV. At the weak scale as is well known, 
the soft terms are no longer universal due to the effects of the renormalization group (RG) 
running. The presence of the RH neutrinos of eq. (2.1) at an intermediate scale contribute 
to the running and generate flavor violating entries in the left-handed slepton mass matrix 
at the weak scale [31]. At the leading order these terms can be estimated to be: 

where Mx represents the GUT scale and Mr u , the scale of the k th RH neutrino, tuq and Aq 
stand for the usual universal soft mass and trilinear terms at the high scale. Y u , the Dirac 
neutrino Yukawa couplings are free parameters in the Type I seesaw mechanism which 
cannot be completely determined even after including the complete data on the neutrino 
mass matrix [32]. 

SO(10) models with their matter representations being 16- dimensional provide a nat- 
ural setting for the seesaw mechanisms. Furthermore, they provide information about the 
neutrino Yukawa couplings. For example, it is known that as long as we restrict to renor- 
malisable SO(10) models, at least one of the neutrino Yukawa couplings should be as large 
as the top Yukawa coupling [25]. Thus with suitable assumptions for the (left-handed) 
mixing of the Dirac Yukawa Neutrino mass matrix, one can make predictions for the flavor 
violation generated at the weak scale from eq. (2.2). Two extreme scenarios for mixing are 
typically considered to be present in Y u [25, 26, 33] : 

Y„ = Y u (CKM Case) 

Y u = Y^ iag Upmns (PMNS Case), (2.3) 

where Y u = VcKMY« iag V^ KM . Both these scenarios can be motivated from concrete 
models of fermion masses within the SO(10) framework [25, 26]. The flavor violating off- 
diagonal entries at the weak scale, eq. (2.2), are then completely determined by assuming 
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Generations 


PMNS 


CKM 


A12 




Y t 2 V td V ts 


A 23 


Y t 2 U^U r3 


Y t 2 V tb V ts 


A 3 i 


Y t 2 U e3 U r3 


Y t 2 V td V th 



Table 1. The dominant combination of neutrino Yukawa couplings which enter eq. (2.2) in CKM 
and PMNS mixing case. 



Yj, as in eq. (2.3). The dominant combinations of Yukawa couplings which enter the radia- 
tive generation of {A^-)ll are shown in table (1). Notice that the hierarchical structure 
of Y u dictated by SO(10) determines that the leading contribution corresponds to third 
generation particles running in the loop. Hence, the flavor violating entries (A12) ll re- 
sponsible for fj, — > e7 process and (A\ 3 )ll responsible for r — > e*y depend on U e3 ~ #13 in 
the PMNS case. The branching ratios of the LFV decays can be roughly estimated to be 



BR{k ^Ijuu) G 2 F raL, 

where m susy is a typical SUSY mass and the flavor violation is as usual parameterize by 
the following quantity 

A f 

4 = -4- (2-5) 

Let us now turn our attention to the NUHM1 boundary conditions. 2 At the first sight 
one might expect that such a modification has no significant effect on the LFV amplitudes, 
except for those due to the modifications in the sparticle spectrum. However, it turns out 
that this is not the only modification. The flavor mixing structure of the slepton mass 
matrix can also be strongly affected. The radiatively generated flavor violating entries in 
eq. (2.2) take the following form in NUHM models: 

where mn u is the soft mass of the up-type Higgs at the high scale. In the present work, we 
consider the NUHM1 scenario, i.e. mn u = m H d at the GUT scale. Furthermore, there can 
be a relative sign difference between the universal mass terms for the matter fields (that 
we still call niQ with abuse of notation) and the Higgs mass terms at the GUT scale. This 
can clearly lead to cancellations (for m\ ~ —2 to 2 ,) or enhancements (for m\ > m o)) m 
the magnitude of the flavor violating entries at the weak scale compared to mSUGRA. 



3 Parameter Range and Phenomenological Constraints 

As mentioned earlier, we will consider two sets of boundary conditions for the soft-terms 
in our numerical analysis. While the mSUGRA is characterized by the standard 'four 

2 Lepton flavor violation in NUHM models has been previously studied in [34], where correlations between 
fi — > e + 7 and /i — > e conversion rates have been discussed. 
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and half parameters (mo, M1/2, Aq, tan/3, sgn(/x)), we parametrize the NUHMl case by 
rriH u = fnH d = "io ~~ Amff. Considering the present and future LHC accessible regions as 
well as the reach of future flavor physics experiments, we scan the soft parameter space in 
the following ranges: 



m G [0, 5] TeV 

f for mSUGRA 

Am H G < 

(JO, 5] for NUHMl 

m 1/2 G [0.1, 2] TeV 

A G [— 3mo, +3m ] 

sgn( M )G{-,+} (3.1) 

Note that we use the convention in which vr? H = sgn(m# u ) |m,#J 2 . For this range of the 
parameter space the first two generations squarks have masses up to 2 ~ 7 TeV and 
the first two generations sleptons up to mj* — 5 TeV. We include in our scan such spectra 
beyond the reach of direct SUSY searches at the LHC, in order to check the capability of 
the flavor violating observables in constraining the parameter space. 

The numerical analysis is carried out using the SUSEFLAV package [35]. It evaluates 
2-loop MSSM RGEs with full 3x3 flavor mixing effects and also incorporates one-loop 
SUSY threshold corrections in all the MSSM parameters. It checks for consistent Radiative 
Electroweak Symmetry Breaking (REWSB) by minimizing the one- loop corrected effective 
superpotential. The program incorporates the effect of RH neutrinos on MSSM RGEs and 
calculates the branching ratios of various LFV processes induced by such RGE effects. The 
program also calculates BR(6 — > s~f) in the minimal flavor violation assumptions. We also 
calculate the BR(B S -> fi + fj,~) using ISABMM subroutine of ISA JET [36]. The light Higgs 
mass is computed using the full two loop corrections of [37-40]. First, we collect the points 
which (a) successfully give REWSB, (b) have no any tachyonic sfermions at the weak scale 
and (c) have the lightest neutralino as Lightest Supersymmetric Particle (LSP). Then we 
calculate all the LFV observables, BR(6 — > 57) using the SUSY spectrum evaluated for 
each point. Finally, we impose the following experimental constraints on the data points 
we collected. 



121.5 GeV <m h < 129.5 GeV 

m^± (lightest Chargino mass) > 103.5 GeV [41] 

BR(B S -»• (J, + fM~) < 4.5 x 10- 9 [42] 

2.85 x 10 -4 < BR(6 -> sj) < 4.24 x 10~ 4 (2<r)[43]. (3.2) 

In comparing our predictions for with the experimental range of eq.(l.l), we take into 
account 3 GeV of theoretical uncertainty (for a recent discussion see [44]). We have not 
considered the Supersymmetric solution to {g — 2)„ discrepancy in the present work. 

In our study, we assume normal hierarchy in the light neutrino mass spectrum and set 

m Ul = 0.001 eV, m V2 = J Am% ol + ml and m vz = J Ara^ tm + m* (3.3) 
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LFV Process Present bound Near future sensitivity 

of ongoing experiments 

BR(/z 67) 
BR(^ — > eee) 
CR(/i -> e in Ti) 
BR(r -> 67) 
BR(r ->■ eee) 
BR(r -»■ ^7) 
BR(r — )■ /i/U/i) 



Table 2. Present bounds and expected sensitivities on LFV processes. 

where AmL and Am^ are the solar and atmospheric squared mass differences respec- 
tively and we use the central values obtained from recent global fits on neutrino data [48]: 

Am 2 ol = 7.62 x 1(T 5 eV 2 and Am 2 tm = 2.53 x 1(T 3 eV 2 (3.4) 

Regarding the mixing angles, we take the most conservative value for the reactor mixing 
angle and set \U e s\ = 0.11, that corresponds to the lower limit of the 3a range given by 
RENO [10] as well as by the global fits [48]. The remaining two angles of the PMNS matrix 
are set to their global fit values #12 = 33.5° and #23 = 45°. The masses of heavy neutrinos 
which we use in our analysis are: 

M Rl = 10 6 GeV, M R2 = 10 9 GeV, and M Ra = 10 14 GeV. (3.5) 

The recent results from the MEG collaboration [8] have improved the existing bound 
on BR(/i — > &y) by one order of magnitude. The present limits on different LFV observables 
are summarized in table 2. In the following subsections, we discuss in detail the results of 
numerical analysis carried out in the PMNS and CKM cases with mSUGRA and NUHM1 
boundary conditions and present a quantitative comparisons between them. 

4 LFV in SUSY SO(10) 

Let us start considering the PMNS case, where there is a direct link between LFV processes 
and neutrino parameters. Besides the CP violating phases, 6\% was the only unknown 
parameter in the leptonic mixing sector till some time ago. There have been various 
theoretical models based on the idea that #13 could be very small, close to zero. In such 
a case, the resulting fi — > e + 7 and r — > e + 7 rates could have been suppressed [49, 50]. 
Recent experiments prove the contrary. Both the experiments measuring #13 are in good 
agreement with each other and indicate a sizeable value of #13 . The recent global fit analysis 
[48] also leads to similar value of 613. All these results indicate that the smallest value of 
\U e s\ is 0.11 at 3d. The resulting rates of BR(/x — > e^) are significantly enhanced for such 
a value of \U e %\. As we are going to see, the updated MEG limit together with a large #13 
puts significant constraints on mSUGRA for moderate as well as large value of tan j3. 3 

3 The interplay between large #13 and BR(/i — > ej) in the context of discrete flavor groups have been 
recently discussed in [51]. 



2.4 x 10" 12 [8] 
1.0 x 10~ 12 [46] 
4.3 x 10~ 12 [47] 
3.0 x 10" 8 [43] 
3.0 x 10- 8 [43] 

4.5 x 10" 8 [43] 
2.0 x 10" 8 43 



10~ 13 [45] 



10~ 8 [45] 
3 x 10~ 9 [45] 
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500 1000 1500 2000 U 1000 2000 3000 4000 5000 

M v 2 [GeV] M [GeV] 

Figure 1. The figure in the left panel shows the BR(/i — > ej) obtained by scanning the mSUGRA 
(in red color) and NUHM (in green color) parameters in the ranges given in eq. (3.1) and for fixed 
tan (3 = 10 and U e 3 =0.11 (the lowest value allowed at 3u by recent RENO observation) and satisfy 
all the constraints in eq. (3.2). Different horizontal lines correspond to present and future bounds 
on BR(/j — > ej). The figure in the right panel shows the allowed space in the mo — m^2 plane 
which satisfy the current MEG bound and eq. (3.2). The region below the red line is excluded by 
the direct searches for SUSY at the LHC [7]. 




M v 2 [GeV] M [GeV] 
Figure 2. The same as figure 1 for tan /3 = 40. 



In figures 1 and 2 we present the constraints from BR(^ — > ej) on mSUGRA and 
NUHM1 parameter space for tan (5 = 10 and 40 respectively. As can be seen, while only 
small part of the paramater space survives for tan /3 = 10 in mSUGRA, it is completely 
ruled out for tan/3 = 40. The allowed regions for low tan/3 require very heavy spectra, 
i.e. mo > 4 TeV for small Mm or My% > 2 TeV for small mo- What is surprising is 
that the constraint on the NUHM1 parameter space is not as weak as one might expect 
form eq. (2.6). As we can see from the figures even in the presence of partial cancellations, 
most of the NUHM1 parameter space is going to be explored by MEG. If one removes the 
light Higgs mass constraint, points with stronger cancellations would be allowed, even with 
[i — > e7 rates below the MEG sensitivity. Thus points compatible with the Higgs mass 
bound, eq. (3.2), do not allow strong cancellations in the flavor violating entry in eq. (2.6). 
For the large tan/3 case, the fi — > e^f constraint is so strong that only few points with 
Mi j2 > 800 GeV are allowed. In the section 4.1, we will discuss in more detail about the 
impact of the constraint on nih in mSUGRA and NUHM1. 

In the context of the updated MEG limit on BR(/i — > ej), it is now worthwhile to see 
what is the situation with the small mixing CKM case. Here we compare the CKM case 
and the PMNS case with mSUGRA boundary conditions. As above, red points correspond 
to the PMNS case while we use the blue color for CKM case. The CKM case has highly 




Mi/2 [GeV] M [GeV] 



Figure 3. The figure in the left panel shows the BR(/i — > &y) obtained by scanning the mSUGRA 
parameters in the ranges given in eq. (3.1) and for fixed tan /? = 10 and U e s = 0.11 (the lowest value 
allowed at 3<r by recent RENO observation). The red (blue) colored points correspond to PMNS 
(CKM) case. Different horizontal lines correspond to present and future bounds on BR(/j — > ej). 
The figure in the right panel shows the allowed space in the mo — mjrt plane which satisfy the 
current MEG bound. The region below the red line is excluded by the current LHC searches [7]. 
Both the plots satisfy all the constraints in eq. (3.2). 



suppressed branching fractions due to the smallness of CKM angles (see table (1)) as has 
been detailed in [26]. Though there has been no strong improvements in the experimental 
sensitivity compared to the analyses of [26], we update the result with the light Higgs mass 
constraint. In figure 3 we show the results for tan f3 = 10. As we can see, some part of the 
parameter space of the CKM case can be probed by the proposed Project-X experiment 4 
for fx — )■ ery. At present the main constraint to this scenario is simply provided by the rrih 
range of eq. (1.1), that excludes the regions with lighter SUSY spectra: niQ < 2 TeV for 
small Mi/2, Mi/% ^5 1 TeV for small mo, as we can see from the right panel of the figure. 
We can also notice that the LHC limits on the mSUGRA parameter space has already 
started to constrain regions of the parameter space otherwise allowed by the bounds in 
eq. (3.2). 

Let us now turn our attention to other observables like [i — > eee, fi — > e conversion 
in nuclei and r — > /17, which is independent of #13. In figures 4, 5 and 6, we show the 
predicted rates for r — > /ry, fj, — > eee and fj, — > e conversion in the Titanium nucleus versus 
the BR(/z —> ej) (that is at present the most constraining LFV observable), for the PMNS 
case in mSUGRA (red points) and in NUHM1 (green points) as well as for the CKM case 
(blue points). 

As can be seen from figure 4, in the PMNS case, the present MEG limit on BR(/i — > erf) 
implies BR(r — > ^7) < 10 -12 , beyond the reach of the proposed experiments. This is a 
direct consequence of the large value of #13 measured by Daya Bay and RENO. In fact, 
from eq. (2.4,2.5) and table 1, we have: 



BR(r -» fjrf) _ \U t3 U„3\ 2 



x BR(r -»• iivv) ps 0(1). (4.1) 



In the CKM case (blue points), the small mixing angle and the bound are such that 
BR(r -> fjL-y) < 1(T 10 . Thus, the scenarios discussed here allow possible signals of LFV in 



4 In appendix A we present a brief summary of all the future experimental facilities related to the flavor 
violating observables discussed in the text. 
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/x — e transition only and evidence for BR(r 
disfavor them. 



^7) at future experiments would strongly 




Figure 4. BR(r -» ^7) versus BR(/i e7) for the PMNS case in mSUGRA (red) and NUHM 
(green), and for the CKM case (blue) for tan/3=10 (left), for tan/3=40 (right). The different 
horizontal and vertical lines correspond to present and future limits on the respective observables. 




Figure 5. BR(^ -> eee) versus BR(/i -> e^) for the PMNS case in mSUGRA (red) and NUHM 
(green), and for the CKM case (blue) for tan/3=10 (left), for tan/3=40 (right). The different 
horizontal and vertical lines correspond to present and future limits on the respective observables. 

Let us now consider the decay fi — > eee. It is known (see e.g. [52]) that in SUSY (with 
conserved R-parity) the dominant contribution to this process arises from the same dipole 
operator responsible for fi — > ej, hence the correlation of the two processes is striking: 

BR(/x ->■ eee) ~ a em x BR(/u — > ej). (4.2) 

Such prediction is consistent with our results shown in figure 5 for tan/3 = 10 and 40. 
The present bound on /i — > eee comes from the SINDRUM experiment at PSI. At present, 
MEG sets a stronger bound, as expected. However the future sensitivity reach of MUSIC 
and Project-X experiments will be able to go beyond the reach of MEG, testing most of 
the NUHM parameter space of our scan. 

Let us now discuss [i — > e conversion in Nuclei, that will represent one of the most 
important probes of LFV in the future. Project X and J-PARC are envisaging facilities 
where conversion on various Nuclei can be studied (for a review, please see [53]). In the 
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Figure 6. CR(^ -> e in Ti) versus BR(^i -> e 7 ) for the PMNS case in mSUGRA (red) and NUHM 
(green), and for the CKM case (blue) for tan/3=10 (left), for tan/3=40 (right). The different 
horizontal and vertical lines correspond to present and future limits on the respective observables. 

present work, we have computed the [i — > e conversion rate in Titanium. The conversion 
rate is as usual normalized by the capture rate of the muon by the nucleus. 

In figure 6, we present our results for CR(/i — > e in Ti) with respect to /i — > e + 7. 
We find that there is a significant spread in the parameter space. This spread is due to 
the existence of cancellations between the penguin contributions at low M±/ 2 an d in low 
tan/3 regions, which has been noted earlier in the literature [33]. Still, we can see that the 
future experiments will be able to test most of the PMNS parameter space and will start 
to constrain the small-mixing scenario (CKM case). 

The above plots have been obtained for the Titanium Nuclei. As in the case of (i — > eee, 
it has been noted that the dominant contributions are from the dipole operators. In such 
a limit, where only the dipole operators contribute, one can easily estimate the conversion 
rate for the new nuclei, X, by knowing its effective charge, Z e jf, the form factor F(q) and 
the atomic number Z and multiplying the conversion rates presented in the above plots by 
the ratio: 



V 4 

Z eff 


\F(q)fZ_ 


X 


Z eff 


F(q)\ 2 Z_ 


Ti 



(4.3) 



The Form factors and Z e ff for various Nuclei can be found in [54]. 

4.1 Combined effect of the Higgs mass and BR(/i — > ej) bounds 

Given the strong constraints from both BR(/j — > ej) and the light Higgs mass one would 
wonder how much of the total parameter space from eq. (3.1) survives in the PMNS case. 
In figure 7 we plot the lightest Higgs mass as a function of mo, M1/2, A$ and tan/3 in 
the three left (right) panels for the mSUGRA (NUHM1) case. In particular, the left panel 
of the third row shows the asymmetric regions in Aq required by the Higgs mass range. 5 
Flavor constraints instead are approximately symmetric in Aq. The last row in the left 
panel shows the constraint on tan/3. As we can see in the mSUGRA case both low tan/3 
(< 5) and high tan/3 (> 20) are ruled out. The constraints in the low tan/3 are due to rrih 
whereas those at high tan j3 are from BR(/i — > ej). 

5 We remind that the gluino-driven radiative effects gives a negative contribution to the stop A-term At, 
so that, if Aq > only very large values could still provide a sizeable \A t \ at low energy. 
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5 1° 15 20 10 20 30 40 50 60 

taqtf tan^ 



Figure 7. Here we present combined regions of parameter space allowed by BR(^ — > ej) and the 
light Higgs mass (m h ), eq. (3.2), on the PMNS case in mSUGRA and NUHM1. 

In mSUGRA/CMSSM models, a light Higgs mass around ~ 125 GeV requires either 
very heavy stops (t± 2) ~ 4 TeV or a large stop mixing triggered by a large low-energy 
value of the stop A-term, \At\. In either of these cases, one can easily convince oneself, 
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using eq.(2.2) that the flavor violating parameter y*i^j) LL * s n °t suppressed. Thus a light 
Higgs mass ~ 125 GeV does not necessarily mean a suppressed flavor violating entry in 
spite of the largeness of stops (£1,2) or ^4-terms required. In fact flavor violation constraints 
are still very strong. 

In NUHM1 case this correlation is somewhat lost due to partial cancellation in the 
flavor violating entry 6 . All values of tan/3 are now allowed and is slightly more sym- 
metric compared to the mSUGRA case. The surprising thing is that imposing the light 
Higgs mass constraint restricts the parameter space to be within the reach of MUSIC and 
Project-X proposals. 

5 Summary and Outlook 

The discovery of a Higgs-like boson at the LHC with a mass close to 125 GeV is one of 
the most significant achievements in particle physics of all times. In the present work, we 
assumed that the particle seen at the LHC is the lightest neutral Higgs scalar of the MSSM. 
We then studied the implications of the observed mass range on SUSY seesaw models along 
with recent improvements in BR(/i — > e + 7) and measurements of #13. We looked at Type 
I seesaw model assuming SO(10) relations between the top Yukawa and the heaviest RH 
neutrino Dirac mass. We assuming two extreme cases of mixing to be present in the Dirac 
Yukawa matrix (a) small CKM like mixing and (b) PMNS like mixing. 

We find that there is a strong complementarity in the PMNS case between the light 
higgs mass constraint and BR(/i — > e + 7). The lower tan f3 regions are strongly constrained 
by the recent measurement of the light Higgs mass and high tan (5 regions by present limit 
of MEG. What is surprising is that tan/3 > 20 is already ruled out. 

The NUHM1 boundary conditions is one more interesting framework where the in- 
terplay between the light Higgs mass constraint and the LFV constraints comes to play. 
To relax the LFV constraints one would need strong cancellations in the flavor violating 
entry, however regions with large cancellations are not favored by a light Higgs mass of 
around 125 GeV. Partial cancellations are however allowed which put these regions within 
the reach of MEG (Project X) for fj, — > e^y (fj, — > eee). 
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A Description of Future Experiments and Prospects: Circa 2020 

PRISM [57, 59] or Phase Rotated Intense Slow Muon source is an upcoming facility at 
J-PARC. It will accelerates muon beam using a magnetic field inside a muon storage ring 

6 In NUHM1 case the cancellations are constrained by the parameter choice of eq.(3.1). 
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LFV process Experiment Future limits Year (expected) 

BR(/x &y) MEG [8] 0(lQ- y6 ) ~ 2013 

Project X [55] O(10" 15 ) > 2021 

BR(/x -> eee) Mu3e [56] 0(1O~ 15 ) ~ 2017 

O(10- 16 ) > 2017 

MUSIC [57] 0(1O~ 16 ) ~ 2017 

Project X [55] 0(1O- 17 ) > 2021 

CR(fi -> e) COMET [57] O(10~ 17 ) ~ 2017 

Mu2e [58] O(10- 17 ) ~ 2020 

PRISM/PRIME [57, 59] ©(lO^ 18 ) ~ 2020 

Project X [55] O(10~ 19 ) > 2021 

BR(r -»• nj) Belle II [60] e>(10~ 8 ) > 2020 

BR(r -> /Ltju/i) Belle II [60] ©(lO^ 10 ) > 2020 

BR(r -> e7) Super B [45] O(10" 9 ) > 2020 

BR(r -»• /i 7 ) Super B [45] e>(10~ 9 ) > 2020 

BR(r -»• /i/i/x) Super B [45] ©(lO" 10 ) > 2020 



Table 3. Future sensitivities of next-generation experiments. 



and will deliver 10 12 ///second at an energy of few 10s of MeV. PRIME [57, 59] is a /i — e 
conversion detector for PRISM designed to stop clumps of \i in a thin foil, having an 
energy of about 20 MeV. Because of high monochromaticity, large pulse rate and very 
low background level it would reach a sensitivity of detecting branching ratios of (fi — > e) 
conversion in nuclei around 10 -18 within few years of running. The PRIME will improve 
the sensitivities to fi — e conversions by two orders of magnitude compared to the next- 
generation experiment COMET (Coherent Muon to Electron Transition) experiment. 

Project X [55], proposed at Fermilab, is a next-generation experiment which has the 
potential to deliver very high power unprecedentedly intense //-beams for precise measure- 
ments of the rare muon decays. Project X will use 0.5 MW beam of \x accelerated at 3 
GeV which will be generated from high-power primary proton beam. With such an intense 
continuous muon beam, Project X can look for \i — > ej with sensitivity of 0(1O -15 ) which 
is two order of magnitude improvement over the future reach of MEG. However, sensitivity 
below 10" 15 appears beyond the reach of Project X unless innovative ideas regarding the 
detectors emerge. Like fj, — > ej, searches for fj, — y 3e also need continuous muon beam. 
Experiments at Project X will improve the sensitivity to [i — > 3e decays by at least three 
to four orders of magnitude, assuming significant efforts to develop new detector technolo- 
gies. The Project X proposals for studying the fi — e conversion in nuclei are particularly 
interesting. Two distinct scenarios has been proposed. If the next-generation round of ex- 
periments (COMET and Mu2e) observe any signal of \i — e conversion, then the available 
muon beams at Project X would allow further precision studies with several fi — e events 
in different nuclei. If no signal of fi — e will be found in next-generation experiments, 
Project X could reach the sensitivities 0(1O~ 19 ) or beyond given some improvements in 
beam technologies. To achieve these goals, a proposal to use muon storage ring installed 
in the muon beam line is under consideration. 

The Belle II [60] is an upgrade of the existing Belle detector at the KEK B-factory in 
Japan. It will make use of the upgraded Super KEKB accelerator. The Belle II detector 
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is expected to collect 40 times more luminosity (8 x 10 35 cm -2 s -1 ) than the previous 
generation Belle detector and by 2022 it would collect an integrated luminosity of ~ 50 ab _1 . 
The sensitivity for the r two body decays is expected to improved by an order of magnitude. 

Super B factory [45] is a proposed high luminosity electron-positron collider. With an 
integrated luminosity of 75 ab _1 , Super B would be able to explore a significant portion of 
parameter space of new physics scenarios by searching for LFV in r decays. It is going to 
improve the sensitivities to different channels of r decays by at least an order of magnitude 
or two. The future sensitivities of all these experiments are summarized in the table 3. 
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